Death of linear response and field-induced dispersion in subdiffusion 
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We discuss the response of continuous time random walks to an oscillating external field within 
the generalized master equation approach. We concentrate on the time dependence of the two first 
moments of the walker's displacements. We show that for power law waiting time distributions 
with < a < 1 corresponding to a semi-Markovian situation showing nonstationarity the mean 
particle position tends to a constant, and the response to the external perturbation dies out. On 
the other hand, the oscillating field leads to a new additional contribution to the dispersion of the 
particle position, proportional to the square of its amplitude and growing with time. These new 
effects, amenable to experimental observation, result directly from the non- stationary property of 
the system. 

PACS numbers: 02.50.-r, 05.40.Fb 



Continuous time random walks (CTRWs) with on- 
site waiting time distributions being power laws lacking 
the first moment have been shown to provide a power- 
ful tool to describe systems which display subdiffusion. 
These subdiffusive CTRWs are non-Markovian (semi- 
Markov) processes characterized by nonstationarity (ag- 
ing). Examples are dispersive charge transport in dis- 
ordered semiconductors, contaminants transport by un- 
derground water, motion of proteins through cell mem- 
branes and many others (see e.g. |X, 3\ for reviews 
and popular accounts) . In the absence of time-dependent 
external perturbations CTRW is a process subordinated 
to simple random walk, thus leading to the description 
within a framework of fractional Fokker-Planck equa- 
tions P, Q) S IE 111- The response to time-dependent 
fields, which has been less explored, is a much more 
delicate problem, due to nonstationary nature of such 
CTRW leading to a variety of effects related to aging 
0, 0, 0i ^1 0|' Here we investigate the response of 
a particle which performs a non-stationary CTRW to 
an external oscillating field by calculation the first two 
moments of its displacement. New effects are observed 
which are absent in the Markovian case. This should be 
also amenable to experimental observation. 

Let us first describe the model adopted throughout this 
work. In a biased decoupled CTRW a particle arriving 
to a site i at time t' stays there for a sojourn time t, 
distributed with the probability density function ip{t) in- 
dependent of external perturbation. This waiting time 
PDF is considered to follow a power-law ip{t) oc TQt~^~°' 
where tq gives the typical timescale for a jump. Leaving 
a site it makes a random step of length a in either di- 
rection. This step is assumed to be instantaneous on the 
timescale of typical waiting times and changes of exter- 
nal parameters; the direction of this step can be biased 
by the time-dependent external force f{t). In what fol- 
lows we turn to dimensionless units and measure length 
in units of a and time in units of tq. The probabilities of 



going to the right iyi-i,i(i) (from site i — 1 to site i) and 
to the left (from site i -I- 1 to site i) are assumed to be 



it) 
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Our description of response of CTRW to external 
fields is based on generalized master equation (GME) 
approach. Several different ways to derive the corre- 
sponding equations are known in the literature (e.g. 
13i,|13i,il^,15, 16J), the one especially fitted to describ- 
ing response to time-dependent fields is given in [T^ . For 
the sake of completeness, we give here the sketch of this 
derivation. The GME follows from two balance condi- 
tions, the probability conservation in a given state and 
under transitions between different states. 

The probability balance for the site k reads 



(2) 



(where the dot denotes the time derivative) with j^{t) 
denoting the gain and loss currents for a site. A particle 
leaving its site k at time t either was in k from the very 
beginning or arrived at k at some Q <t' <t so that 

j-(t) = mpk{o)+ 1 ^{t^t')]+{t')dt' (3) 

Jo 

= mpk{<^) + I v^t - 1') [pk{t') + j^it')] dt\ 

Jo 

where in the second line Eq.(|2Jl was used. The formal 
solution to this equation can be expressed through an 
integro-differential operator 

j^it) = $p,it) = ±j\lit^t')p,it')dt' (4) 

with the memory kernel M{t) given by its Laplace trans- 
form 



M{u) 



i-i}{uy 



(5) 
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The probability conservation for transitions between 
different sites give the relation between the gain current 
in the state k and loss currents at neighboring sites: 



(6) 



Inserting the corresponding expressions into the first bal- 
ance equation gives a GME for Pk{t): 

Pk{t) ^Wk-i,k{t)^Pkit) + Wk+iM{t)^Pkit)-^Pkit)- (7) 

Note that the integro-differential operator $ does not 
commute with the function of time Wij{t). 

Using now Eq.(^ for the transition probabilities and 
passing to the continuum limit we get a generalized 
Fokker-Planck equation 



(8) 

For Markovian random walk process with exponential 
waiting time distribution (corresponding to a = 1 and 
thus to M{t) = 1) this equation reduces to a usual 
time-dependent Fokker-Planck equation. For power-law 
waiting-time distributions t^^~°' with < a < 1 

one gets M{t) oc t"^^ and the integro-differential op- 
erators on the right-hand side of this equation get pro- 
portional to the operator of fractional Riemann-Liouville 
derivative ^ M{t ~ t')f{t')dt' cx QDl-"f{t). This is 
exactly the case we now concentrate on. 

The CTRWs with < a < 1 are known to show a va- 
riety of phenomena connected with non-stationarity (re- 
lated also to the so-called aging property 0i S Qi 
lio, .21.] ) . One of the manifestations of aging is the decay 
of response of the system to an alternating or pulsed field 
in course of the time (i.e. when its age grows), see Ref. 
1^ . Here we consider the response of the system to a sinu- 
soidal external field f{t). We start from Eq.ljSj) and con- 
sider the moments oi p{x, t), ninit) — x^p{x, t)dx of 
the probability distribution of particle's positions. These 
moments can be easily obtained by multiplying both sides 
of Eq.® by a;" and integration. Assuming the system to 
be infinite and spatially homogeneous we get by partial 
integration 



x^^\/p{x,t')dx — —mrin-iit') 



(9) 



and 



'Ap(x, t')dx = n{n - l)TO„_2(i') (10) 



(for n > 2). Thus, general equations for the moments are 
given by 

TliTl — 1 ) 

mn{t) = nfifit)^mn-i{t) + ^$m„_2(t). (11) 



To be able to use these equations in the whole range 
oi < m < oo one can formally put mo = 1 and 
TO_i(t) = 0. The equations for the first moment (mean 
displacement) and the second (dispersion) read: 



and 



m,it)=^ifmi = ^ifit)M{t) 



m2{t) = 2^/(t)$TOi(t) + M{t). 



(12) 



(13) 



Note that for semi-Markovian cases with < a < 1 
M{t) cx t"^^ is the decaying function of time. Therefore 
the response mi {t) to the perturbation vanishes in course 
of the time leading to the effect we call "death of linear 
response" in systems showing subdiffusion. 

The second moment, m2{t) = /g m2(t')dt' consists es- 
sentially of two contributions, m2{t) = (7i{t) -l-cr|(t), the 
one depending on the external perturbation (through the 
first moment mi) 



af{t)=2 dt'iif{t')— M[t' ~t")mi{t")dt" (14) 

and of the field-independent purely (sub)diffusive contri- 
bution 



crl{t) = / M{t')dt' 



(15) 



Let us discuss the overall structure of expressions for 
the first and the second moment for the case of a periodic 
force f{t) — fosinujt. The first moment is: 



miit) 



/i/o / dt' s\nujt'M{t')dt' 
Jo 



dt'- 
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-M{t')dt'. (16) 



Turning to the Laplace domain we get 

mi[u) = — \m{u - iu) - M{u + iu) 
2iu L 



(17) 



as it follows from the shift theorem. The asymptotic be- 
havior of the first moment then follows straightforwardly. 
For u ^ (corresponding to i — s- oo) we have 



rhi{u) 



2iu 



M{-iuj) - M{iuj) 



(18) 



being the Laplace transform of a constant mi(oo) = 
—fifolmM{—iLu). The Laplace transform of the field- 
dependent contribution to the second moment (again ob- 
tained by using the shift theorem) reads: 



aliu) 



4w 



|M(M-iLj) M{u-2iLu) — M(u) 



M{u + ioj) M (u) - M{u + 2iuj) | 



(19) 
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To obtain the asymptotic behavior of (T^(m) it is enough 
to note that for power-law waiting time distributions with 
< a < 1, M(u) oc which diverges when u 0. 
Thus, the leading contribution to o'^(m) is 



2 2 M{-iu))M{u) + M{iuj)M(u) 



Au 
M{u) 



(20) 



This means that the asymptotic growth of the field- 
dependent contribution to the dispersion is given by 

The most important feature of this result is that al- 
though the first moment of the distribution stagnates, 
the field-dependent contribution to the dispersion contin- 
ues growing, a manifestation of a new effect, specific for 
non-stationary CTRWs, namely the field- induced disper- 
sion. This growing contribution is absent for the Marko- 
vian and asymptotically Markovian processes (a = 1) 
only due to the fact that the corresponding prefactor 
Re(iw)~" vanishes. 

To get an impression on the overall behavior of the first 
and the second moments let us consider the special case 
a = 1/2. In our numerical example we set w = 1. In this 
case 



d r* 



(21) 



= r(l/2)Ai/osini o^y'l = M/oi"'^'sini. 

Note that the fractional derivative of the constant does 
not vanish. From this expression we see that the sus- 
ceptibility of the system to an external force decays in 
course of the time as i^^/^; namely its linear response 
dies out. Integrating this expression over time we get the 
analytical form for mi{t) 



(22) 



involving the Fresnel sinus-integral S{x) = 
(l/-\/27r) dt sin (t) /-s/t. The Fresnel integral tends to 
a constant value for large values of its argument; the 
final value of the first moment is mostly determined by 
the value of the external perturbation at short times, 
when the system was "young" , the result of what was 
called "Freudistic" memory of aging systems in Ref.j^. 
The behavior of mi{t) is shown in Fig. 1. 

The behavior of the field-dependent contribution to the 
second moment is 



ol{t)^^ihl dt'sint'-^ I {f -t")-"^mAt")dt 



"'0 



Integrating this expression by parts twice we get: 



sm{t)^{t)- / ■${t')cos{t')dt' 
Jo 



(23) 



(24) 
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FIG. 1: The mean displacement mi{t) (measured in units of 
fJ^fo) in a CTRW model with a = 1/2 as a function of time 
for uj — 1. Note that for t large the displacement stagnates: 
The linear response of the system to external sinusoidal field 
dies. 



with 



*(0 



^y{t~t')/t' aint'dt' 



(25) 



The numerical evaluation of the integral is shown in Fig. 2 
showing [crj(t)]^/", i.e. the square of the corresponding 
expression, Ea. (|24|l . The leading contribution t" to the 
overall behavior corresponds to the overall linear growth 
of [af{t)]^/°'. Interesting is the subleading asymmetric 
oscillatory behavior overplayed on this overall growth. 
The amplitude of these oscillations decays on a linear 
plot, but seems to stay constant in the one presenting 
the square of the function. This means that the decay of 
the subleading term follows essentially t^^^^. The field- 
independent subdiffusion contribution cr|(t) grows as t^/^ 
according to Ea. l(TB|) . 

Let us summarize our findings. We discussed the be- 
havior of a particle performing continuous-time random 
walks with a power-law distribution of waiting times lack- 
ing the first moment {tpit) cx t^^^" with < a < 1) un- 
der the influence of oscillating external field. Using the 
approach based on the generalized master equation we 
derive equations for the first two moments of the displace- 
ment. The first moment of the displacement stagnates, 
an effect we term "death of linear response" . The sec- 
ond moment, on the contrary, grows as and contains, 
in addition to the normal (sub-) diffusion contribution, a 
field-induced contribution, proportional to the square of 
the external field. This new effect which shows up in 
non-stationary CTRW is absent in the Markovian case 
(a = 1) since the corresponding prefactor vanishes. 
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FIG. 2: Shown is the square of ai{t) (measured in units of 
fa) as a function of t for a = 1/2 and ui — 1. 
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